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ABSTRACT 

We show by means of a high-resolution N-body simulation how the mass assembly histories of 
galaxy-size cold dark matter (CDM) halos depend on environment. Halos in high density environments 
form earlier and a higher fraction of their mass is assembled in major mergers, compared to low density 
environments. The distribution of the present day specific mass aggregation rate is strongly dependent 
on environment. While in low density environments only ~ 20% of the halos are not accreting mass at 
the present epoch, this fraction rises to ~ 80% at high densities. At z = 1 the median of the specific 
aggregation rate is ~ 4 times larger than at z — and almost independent on environment. All 
the dependences on environment found here are critically enhanced by local processes associated to 
subhalos because the fraction of subhalos increases as the environment gets denser. The distribution of 
the halo specific mass aggregation rate as well as its dependence on environment resemble the relations 
for the specific star formation rate distribution of galaxies. An analogue of the morphology-density 
relation is also present at the level of CDM halos, being driven by the halo major merging history. 
Nevertheless, baryonic processes are necessary in order to explain further details and the evolution of 
the star formation rate-, color- and morphology-environment relations. 

Subject headings: cosmology:dark matter — galaxies:formation — galaxies: halos — methods:iV— body 
simulations 



1. INTRODUCTION 

The build-up of galactic dark matter halos is a crucial 
ingredient of galaxy formation in the context of the pop- 
ular hierarchical A Cold Dark Matter (ACDM) scenario. 
Many of the present-day galaxy properties are expected 
to be tightly related to their halo mass assembly pro- 
cess. A simple way to characterize this process is by the 
halo mass-aggregation history (MAH). The MAH im- 
plies both the mass growth by violent major mergers and 
by quiescent accretion. The influence of the MAHs on the 
halo and galaxy properties re mains imprinted for exam- 
ple on the halo concentration llA vila-Reese et al. Nl998t 
iFirmani fc Avila-Reese l2000tlWechsler et al. I2002|) . the 
gas infal l and star formation rates in disk galax- 
ies (e .g.. lAvila-Reese fc Firmani 1 120001 Ivan den Bosch I 



2002 : IMuralietalJ 2002) and the maj or merging rate and 
there f ore, the morphology of galaxies (iK^fhnMin eiLal 
1993 iBaugh et al I ll9M iSomeryille fc Primackl 1199' 
" Cole et al. _ l l2000t ISprineel et all 120011 iGranato et al 



present-day mass, starting from a primordial Gaussian 
density fluctuation field characterized by a power spec- 
trum. On one hand, the results have shown that, due 
to the stochastic nature of the primordial density field, 
the MAHs of halos of a given mass at z — span a 
wide range of "tracks" , influencing galaxy features such 
as the scatter of the Tully-Fisher, halo concentration- 
mass and color-magnitude relations ( Eiscnstcin_fclj0^b_ 
1996; jAy^kvP^ge^^t^d^lP^UFm 
| 2000tlBullock et al. 1120011 lEke et al.ll200lUBerlind et al 
2005). On the other hand, the MAHs depend on mass: 
less massive halos tend to build up a given fraction of 
their present-day mass on average earlier than the mas- 
sive ones. This dependence is at the basis of the halo 
conce ntration-mass relation i j Na vjirro. F^nk^ Whi^gl 
T9971 IFirmani fc Avila,-Reese I 1200(1 . lEke et al 



2003b). 



van den Bosch H2002t IWechsler et al. 1120021 iZhao et al I 



200lHSteinmetz fc Navarroll2002HMaller et al. 112005^ 

The continuous mass growth of isolated halos is a 
generic proce s s in the hierarchical CDM scenario (e.g., 
[ Gunn I 119821: iKauffmann et al. I 11993 lLacev fc Colel 
1993). Semi-analytic models and numerical simulations 
have been used to predict the MAHs of halos of a given 
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To study the general behavior of the MAHs as a func- 
tion of mass and environment, we introduce the average 

MAH (AMAH) as M(a) = < ^ > , where a is the scale 
factor and < M(a) > is the average mass of the most 
massive progenitors (MMP) of present-day halos of mass 
Mq at the epoch a. It was found that the smooth AMAHs 
can be approximated by simple (universal) functions, 
where (1) the main function parameter is related to some 
typical formation epoch of the halo, and (2) this epoch 
depends on Mq dBullock et al. 1 120011 Ivan den Bosch I 
12001 IWechsler et al. 1 120021) . On the basis of N-body 
cosmological simulations, it was also suggested that all 
the MAHs, independent of Mq, present an early phase of 
fast mass aggregation (mainly by major mergers) and a 
late phase of slow aggregation (mainly by smooth mass 
accretion) (|Zhao et al. Il2003al ISalvador-Sole et al.ll2005l: 
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iLi et alj |2005). These authors also show that the depth 
of the potential well of a present-day halo, traced by 
its virial velocity, is set mainly at the end of the fast, 
major-merging driven, growth phase. 

The role that the environment plays on the CDM halo 
assembly has not been explicitely explored in the litera- 
ture. Yet, it is well known that some galaxy properties 
vary as a function of environment. Present-day galaxy 
properties as well as their evolution have been studied by 
some authors by combining cosmological N-body simu- 
lations and semi-analytic modeling of the baryon galaxy 
proces ses, a technique pi oneered by iKauffmann et all 
(1999) and iSoringel et al] {2001). Although in the pa- 
pers that used this technique the galaxy dependencies 
on environment are intrinsically taken into account, the 
impact the MAH dependence on environment has on the 
galaxy properties is not clearly established. 

The main observable dependencies with envi- 
ronment are seen for the morphological mix of 

elliptical and spir al gala xi es (e.g.. [Dressier I [19801 
Postman fc Geller I 119841: iWhitmore et al. I Il993t 
Dommeuez et al. I 120021 iGoto et al. I I2003T). color 
and s pecific star formation rate ( e.g.. iBalogh e t al. I 
19981 12004 iKodama et al. I 120011 iTran et al. I 120011: 
Pimbblet et al. 1120021 iLewis et al. Il2002t iBlantqn et all 
20031 iGomez et al I 12003 IKauffmann et al. I 120041: 
Hogg et al. I 12004 iTanaka et al. I 12004 ICroton et al. 1 
2005). Besides, there is some evidence that these 
depen dencies evolve, becoming weaker at h igher red- 
shifts llTreu et al. 1 120031: IGoto et al. 1 121101 iBell et all 
12004 iSmith et all 120051: [Postman et al. 1 12005|) . As 
mentioned above, the morphology of galaxies as well as 
their colors and star formation rates are certainly related 
to the assembly history of their halos. Astrophysical 
external mechanisms, acting mainly in high density 
environments, are also important. To disentangle the 
role of one from the other and to understand which 
are the drivers of the changes of galaxy properties with 
environment, it is important to explore and quantify how 
the CDM halo assembly history depends on environment. 
N-body cosmological simulations are required for this 
endeavor. 

Bas ed on N-body simulations. iLemson fc Kauffmannl 
(1999) concluded that only the halo mass function varies 
with environment. No significant dependence of any 
other halo property on environment was f o und b y these 
authors. iGottlober. Klvpin fc Kravtsovl (|2001f ) found 
that the major merging rate histories of CDM halos 
vary as a functi o n of environment. More recently, 
lAvila-Reese et all l)2005t) found that some properties of 
halos of similar masses (for example, the mean con- 
centration and spin parameter) actually vary systemat- 
ically between voids and clusters. Also, recently, it was 
shown that halos o f a given mass form stati s tically earlier 
in denser region s ( Shct h fc Torm^ 120041 fHarker et all 
120061 iGao et al J [20051 iWechsler et alJl2005|) . Nonethe- 
less, the question of how the mass assembly history of 
CDM halos of similar masses depend on environment 
has not been yet explored in detail. Furthermore, we 
would like to know to what extent this dependence is 
able to explain the observed dependencies of galaxy mor- 
phology, color, star formation rate, etc. on environment. 
We would also like to know if the bimodality of the star 
formation rate (SFR) and color distributions, found in 



large statistical galaxy samples, could be accounted for, 
at least in part, by the physics of dark matter halos only. 

In this Paper we construct the MAHs of ~ 4700 
halos with present-day masses larger than 10 11 h~ 1 MQ 
from a 50/i -1 Mpc box simulation with high-mass res- 
olution: m p = 7.75 x 10 7 /i _1 M o (§2). We find that 
several quantities that characterize the halo MAH (the 
specific mass aggregation rate and merging mass frac- 
tion histories, formation times, etc.) change significantly 
with environment for halos of similar present-day masses 
(§3). We also explore possible systematical dependences 
of observational-related quantities on intermediate-scale 
density, and discuss to what extent these dependences 
are able to explain the observed galaxy property-density 
relations (§4). Finally, we highlight the main conclusions 
of our work (§5). 

2. THE SIMULATION 

We adopt a flat cosmological model with cosmologi- 
cal constant (ACDM) and the parameters r2 tllj o = 0.3, 
Oa = 0.7, and h = 0.7. The initial matter power spec- 
trum has been calculated using the numerical results of 
direct integration kindly provided by W. Hu, and it was 
normalized to erg — 0.9. The study in this Paper relies on 
the results from a simulatio n of box size 50 /i~ 1 Mpc run 
with the GADGET-II code (|Springe]|l2005j) . With 512 3 
particles and a particle mass of m p = 7.75 x lO 7 /i _1 M0 
we are able to determine MAHs that extend from z — 
to redshifts as high as 6 — 9 for halos more massive 
than lO n /i _1 M (containing > 1300 particles). In 95 
timesteps of Aa = 0.01 (~ 100—160 Myr), halos are iden- 
tified by a new MPI version of the Bound Density Max- 
ima a lgorithm originally introduced by iKlvpin et al. I 
(1999). This algorithm allows to detect isolated or 'par- 
ent' halos (self-bound structures not contained within 
larger ones) as well as subhalos (self-bound structures 
contained within larger ones). From these halo cata- 
logues two merger trees were constructed, one for isolated 
halos alone and one for all halos, including subhalos. At 
z = there are more than 4700 halos more massive than 
10 11 /i- 1 M o . 

A halo at redshift z\ with nj(z\) particles is identified 
as the MMP of a halo at zq < z\ containing rii(zo) par- 
ticles if at least a fraction f m i n = 0.2 of its particles are 
found in the progenitor and the overlap of particles n ov = 
rii(zo) nrij(zi) divided by n max = maxlni(zo),nj(z{)) 
is maximal, fmin is the only free parameter which is 
chosen to allow for major mergers and early rapid mass 
growth. MAHs and merger trees only very weakly de- 
pend on this parameter. For all progenitors which are 
merging in a timestep, at least half of their particles are 
required to be found in the descendant at zq. The algo- 
rithm to construct full merger trees from simulations will 
be described in a forthcoming paper. The construction 
of MAHs merely requires the identification of the MMP 
which is straightforward. 

The intermediate-scale environment associated to a 
given halo is defined as the density contrast of dark mat- 
ter in a sphere of radius R around the halo center of mass, 
5(R) = p(R)/ph g — 1 (with pbg the background density 
n p crit ). The qualitative conclusions are not dependent 
on the choice of R for 2 < i?//i _1 Mpc < 8. Since the 
dependences on environment become weaker for large R, 
we will use here 64 = d(R = 4/i~ 1 Mpc). Note that with 
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Fig. 1. — The upper panel shows the number distribution of 
halos with density in the 3 mass ranges indicated in the plot. The 
lower panel shows how the fraction of halos in the given mass range 
that are subhalos depends on density contrast. The two vertical 
lines indicate our cuts in density and the definition of low density 
(84 < 0) and high density (£4 > 5) environments. 



the 84 criterion, the local environment is smoothed out. 
For example, the value of 84 for most subhalos contained 
within a given parent halo is the same, notwithstanding 
whether the subhalo is in the center or in the periphery 
of its parent halo. We are here interested in exploring 
the general effects of environment on the halo assem- 
bly process. The more local environmental effects are 
related only to sub halos and have been studied in de- 
tail previously (e.g.. iKravtsov et afll2004bt iReed et al. I 
120051: iDe Lucia et all 112004^ In all cases the density 
contrast is defined at z = 0, except for those panels in 
Figures 8 and 9, where results for z = 1 are shown and 
84 is defined correspondingly at this redshift. 

The distribution of the density contrast 84 is shown in 
the upper panel of Fig. ^ for all halos in the three mass 
bins used hereafter. The vertical lines at #4 = 5 and 
84 = illustrate our definitions of 'high' and 'low' den- 
sity environments, respectively. Density contrasts 84 > 5 
correspond to cluster environments, while 84 < corre- 
spond to the outskirts of filaments and to voids. In the 
lower panel we show the fraction of subhalos in subsam- 
ples selected by mass as a function of 84. Whereas in 
the low density environment only around 1% of the ha- 
los with their maximum mass larger than 10 /i M are 
subhalos, this fraction at high density amounts to about 
half of all halos. The overall fraction of halos more mas- 
sive than 1O 11 /i _1 M0 that are subhalos rises from 7% at 
z = 1 to 14% at z = 0. 

We also explored the density measure used by 

iLemson fc Kauffmann I l)1999j) . which excludes the mass 
of the halo itself and the close neighborhood, and uses 
only matter in a spherical shell of inner radius 2/i _1 Mpc 
and outer radius 5/i _1 Mpc to calculate the density con- 
trast 82-5- The results presented in the next section 
virtually do not change when using 82-5 instead of 84. It 
should be noted that, according to hybrid N-body /semi- 
analytical models, the galaxy n umber density is propor - 
tional to the mass density (e.g.. lKauffmann et al.l ll997L 



3. RESULTS 
3.1. Mass aggregation histories 

In Fig. 0we present and compare the dependences of 
the AMAH on both mass and environment (84). The 
AMAHs of isolated halos (left four panels) and all ha- 
los (isolated+subhalos, right four panels) are shown in 
panels A and E for three different mass ranges: 11 < 
logMo/Zi^Mg < 11.5, 11.5 < logM o //j- 1 M < 12 and 
12 < logMo//i _1 M0 < 13. Further, the AMAHs for 
the three mass bins, each one divided in our two ex- 
treme density contrasts are shown. Dashed and dotted 
curves correspond to the AMAHs in high- {84 > 5) and 
low- (84 < 0) density environments, respectively. These 
density-dependent AMAHs are shown with mass increas- 
ing from panel B to D for isolated halos, and from panel 
F to H for all halos. 

For the construction of the AMAHs we adopt the fol- 
lowing definitions. Individual MAHs are normalized to 
their maximal mass Mi^ rnax . For most isolated halos 
the maximum mass is reached at the present day epoch, 
Mi iTnax = Mi(z = 0). If a halo suffers mass loss, we 
keep Mi(a) fixed until its mass further grows. Therefore 
all AMAHs have M(a = 1) = 1. subhalos actually can 
experience substantial mass loss once they fall into their 
parent halo. However, it is reasonable to assume that 
the mass of a galaxy living in a subhalo is proportional 
to the maximum mass reached by the subhalo rather than 
to the current stripped halo mass: gas accretion related 
to the halo mass growth is stopped after the halo falls 
into a host halo; furthermore, halo tidal stripping is not 
expected to substantially affect the mass of the galaxy 
formed at the core of the halo. Therefore the mass of a 
halo is kept fixed to the mass the halo had at the time it 
becomes a subhalo (M max ), except for those halos that 
suffer mergers within the host halo, increasing thereby 
their mass. 

When individual MAHs are followed to high redshift, 
the mass of more and more halos falls below the resolu- 
tion limit. To avoid a bias in the AMAHs due to this 
incompleteness, the AMAHs are calculated only down 
to timesteps when more than 90% of individual MAHs 
can be followed. This is the reason why e.g., for the 
lowest mass bin, the curves already end at M ~ 0.1. 
Since the slope of the mass fu nction becomes steeper 
in lower density environme nts l)Lemson fc Kauffmann I 
1999; iGottlober et al.ir2003|) . their average mass in a 
given mass bin is lower than for halos in the same mass 
bin but in high density environments. To exclude this 
mass effect when comparing different environments the 
high-density environment halos are selected in their mass 
bin such that the average mass is approximately equal to 
the bin average mass of halos in the low-density environ- 
ment. To achieve this we have to exclude a few of the 
most massive halos in each 84 > 5 bin. 

Figure [21 shows the well known dependence of AMAHs 
on mass: low mass halos on average assemble a given 
mass fraction earlier than massive halos. To get an im- 
pression of the scatter in individual MAHs, the 10 and 90 
percentiles of the distribution are shown by dot-dashed 
lines for the middle mass bin (11.5 < logMo//i~ 1 MQ < 
12). This scatter is much bigger than the differences be- 
tween AMAHs corr esponding to the different mass bins, 
a fact first noted by lAvila-Reese et al~l l)1998|) . 
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Fig. 2. — Average MAHs for isolated halos are shown in the left 4 panels whereas AMAHs in the right 4 panels are built from all halos, 
including subhalos. The halo population is divided into 3 mass bins as indicated in the plot. In panels A and E the AMAHs for all halos 
in these mass bins are compared. The 10 and 90 percentiles for the middle mass bin are shown as dot-dashed lines. From the halo samples 
in the three mass bins two subsamples are selected at low and high density, respectively. These are shown with mass increasing from panel 
B to D and from panel F to H. 



The new result presented here is the dependence of the 
average halo MAHs of similar masses on environment. 
The lower is the density contrast #4, the later on average 
halos of similar final mass accumulate their mass. This 
ambicntal dependence is already seen for isolated halos 
(left panels of FigEJl . However, the dependence becomes 
much stronger when subhalos are included, as shown in 
the right four panels of FigEl As we saw in Fig. ^ 
the fraction of halos that are subhalos is a strong func- 
tion of environment. On one hand, most subhalos reach 
their maximum mass at the time they fall into their host 
halo and have truncated MAHs. On the other hand, the 
properties of subhalos do not depend si gnificantly on the 
host h alo properties or its environment l)De Lucia et al.1 
2004a). Thus, the main contribution of the MAHs of 
dark matter halos to the observed environmental depen- 
dence of galaxies should result from the dependence of 
the subhalo fraction on environment. 

3.2. Mass aggregation rates 

The MAHs can be further characterized by the mass 
aggregation rate of the halos. We show in Fig. [3] aver- 
ages of fractional aggregation rates per Gyr normalized 
to the maximum mass of the given halo, M/M max . In 
building the average M/M max , mass loss is not taken 
into account, i.e. all halos with dM/alt < have their 
aggregation rate set to for the reasons outlined above. 
The left upper panel compares M/M max for the same 
mass bins as in Fig. While at high redshift the low 
mass halos are the ones that accumulate the mass faster, 
at low redshift we find that the high mass halos have the 
higher mass aggregation rates. Thus, the maximum of 
M/M max moves slightly to lower redshifts with increas- 



ing mass. However, the differences in the M/M max his- 
tories with mass are not so dramatic as with environment 
(see the rest of the panels of Fig. |3J|. The average ag- 
gregations rates at z sa are much higher in low-density 
environments than in the high-density ones, but these 
differences are reversed at redshifts higher than z ~ 1. 
The former is mostly due to the presence of subhalos 
which can not accrete more mass after falling into their 
host halo (the only way to increase the mass is through 
mergers but mergers are very unlikely inside a larger viri- 
alizcd structure). 

The distribution of the present-day specific mass ag- 
gregation rate, (M/M max )o, consists of two parts: 43% 
of all the halos in the mass range 11.5 < log A/o/Zi^'Mq < 
12.5 have (M /M max )o < ('passive' halos), while the 
rest has (M /M max )o > 0, with a distribution that peaks 
at ~ 0.04Gyr _1 (See Fig. EJ. The distribution of 
(M/M max )o changes only little with mass, in the sense 
that massive halos have a slightly higher rate than the 
less massive halos (upper panel of Fig. 0}. On the other 
hand, the distribution changes dramatically with envi- 
ronment (lower panel of Fig. QJ- In the mass range of 
11.5 < logMo//i _1 M < 12.5, the fraction of 'passive' 
halos is 81% and 22% for the high- (64 > 5) and low- 
(<?4 < 0) density environments, respectively. The distri- 
bution of (M/M max ) for halos with (M/M max ) Q > 
('active' halos) also depends on environment: the me- 
dian of the distribution is at (M /M max ) w 0.02 Gyr -1 
and ss 0.04 Gyr" 1 for the high- and low- density envi- 
ronments, respectively. 

Our results are not complete if we do not estimate the 
accuracy in the measurement of the mass aggregation 
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Fig. 3. — Evolution of the average specific mass aggregation rate 
per Gyr for halos in different mass ranges and different environ- 
ments. All (isolated+sub-) halos are included here. Line styles are 
as in Fig. [2] 
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rates and the fractions of passive/active halos. The for- 
mer measurement is done by the estimate of the halo 
mass at two different output times. Therefore, its ac- 
curacy is related to the accuracy of the measurement of 
the halo mass in the simulation. The main source of 
error can be due to the halo finder (BDM). Since the 
initial seeds are chosen randomly and then moved it- 
eratively to the centers of candidate halos, the result- 
ing radii and masses of halos can in general be differ- 
ent. Differences due to a changing particle configuration 
can only be detected if they are greater than the dif- 
ferences introduced by the random seeds in BDM. We 
rerun the halo finder on the last timestep with different 
sets of initial seed centers for the halos. If Mi and M 2 
are two different mass determinations of the same halo, 
we find for the variance of (Mi — M 2 )/Mi for all halos 
in the mass range 11.5 <log(M//i~ 1 M Q )< 12.5 the value 
0.00225, and hence for the error in (M /M max )o a value 
of 0.0045. Therefore, except for the two lowest bins with 
\og(M/M max ) < -2.3, Fig. 4 is not affected by these 
errors. The fraction of halos with < (M /M max )o < 
0.0045 is 1.4%. However, we recall that the MAHs we 
are using in our analysis are not allowed to decrease 
with time. If we use the actual MAHs instead, we find 
that ~ 5% of the halos in the above mass range have 
-0.0045 < (M/M max ) < 0.0045. We consider that this 
percentage represents the maximum error in the passive 
halo fraction determination. 

Related to the mass error, there is also a systematical 
uncertainty in our estimate of the passive halo fraction. 
As explained and justified above, we used MAHs that 
are kept constant if the actual mass is falling. When 
using the actual mass in the two timesteps z = and 
z = 0.1, the fraction of passive halos is 36% instead of 



Fig. 4. — The distribution of the aggregation rates of dark matter 
halos between redshift z = 0.1 and is shown for three different 
mass ranges in the upper panel and for halos in the mass range 
11.5 < log Aj r n/'i~ 1 M0 < 12.5 at different density contrasts in the 
lower panel. The fraction of halos with no aggregation or mass loss 
is indicated in the figure. The distribution is shown for all halos, 
including subhalos. 

43%. On the one hand, the actual mass does not suffer 
from a temporary wrong assignment of a too high mass 
value in the past, as does the 'non-decreasing' MAHs. 
On the other hand, the actual mass gives not necessarily 
a better estimate of the passive halo fraction, since the 
actual mass is ignoring the possibility that the halo could 
have reached its maximum mass in the past and therefore 
is passive at present. 

Summarizing, our method in general allows us only 
to find a possible range for the fractions of pas- 
sive halos. For all the halos in the mass range 
11.5 <log(M/ft,- 1 M©)< 12.5 this lies between (36 ± 5)% 
and (43 ±5)%. 

3.3. Formation times 

Another way to characterize individual MAHs is by 
their formation time. Here we use the time when 
the MMP in the MAH reaches half of its maximum 
mass. To isolate the ambiental effect from the known 
mass differential effect we use as variable for the for- 
mation time the scaled formation redshift, ll>i/ 2 . This 
quantity "absorbs" the dependence of z 1 / 2 on mass 
(|Lacev fc Cole I Il993j) and it is defined as LO\j2 = 
[S c (z 1/2 )-S c (z )}/[^/a 2 (M /2) - a 2 (M )], where 5 c (z) is 
the critical density threshold for collapse at z and a 2 (M) 
is the linear theory variance of density fluctuations at 
mass M. 

In Fig. [5] the distribution of 0)1/2 f° r isolated halos 
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Fig. 5. — Distribution of formation times in terms of uJi/2- The 
upper panel uses isolated halos only. The lower panel shows the 
distribution of all halos, including subhalos, in the mass range 11 < 
logM o //i- 1 M <13 



(upper panel) and all halos (isolated + subhalos, lower 
panel) in the mass range 11 < log Mo/h Mq < 13 (his- 
tograms with solid lines) are plotted. We further divide 
both samples in the two populations of high-density en- 
vironment halos (#4 > 5, histograms with dashed line), 
and low-density environment halos (#4 < 0, histograms 
with dotted line). There is a shift in the distributions 
(highly significant according to a K-S test), the former 
having systematically higher scaled formation redshifts 
than the latter. Again, the difference becomes more pro- 
nounced when subhalos are included, because these halos 
live in the densest regions which collapse earlier and the 
mass aggregation is stopped as soon as they fall into their 
host halo. The average values of Co 1/2 for the sample in- 
cluding subhalos are 1.13 ±0.6 and 1.74 ±0.9 for the low 
and high-density environments, respectively. This trans- 
lates to formation redshifts of Zf t i ow = 1.14 ± 0.5 and 
z fMg h = 1.67 ±0.8, respectively for M = 1O 12 /i _1 M . 
In cosmic time, this difference corresponds to ~ 1.5 Gyr, 
suggesting that galaxies formed in the high-density en- 
vironment will be redder than those formed in the low- 
density one. 

For reference we also show in Fig. [5] the pre- 
diction of the exten ded Press-Schechter (EPS) theory 
ijLacev fc Cole lll993h using the ACDM power spectrum. 
Compared to the simulation result for isolated halos (see 
solid curve, top panel), there is a shift to later formation 
times in the EPS app roach as was previously noted by 
Ivan den Bosch I $M)2). 



3.4. Major mergers 

Major mergers are believed to play a crucial role 
in shaping galaxies and leading to an environmen- 
tal dependence of galaxy morphology. We study 
here further the question of the halo major-merging 
dependence on environment. E arlier results by 
iGottlober. Klvoin fc Kravtsov I l)2001l) showed that such 
a dependence exists. We define a major merger when 
the relative mass increase (MMP complement), AM/M, 
is larger than 0.2 in a timestcp of Aa = 0.01. Results 
do not change qualitatively by fixing another reasonable 
mass fraction increase. We choose this criterion to be 
able to find mergers at high redshifts, where AM can 
still be determined for small halos but no secondary 
progenitor would be found. However, for subhalos we 
have found that the number of major mergers is over- 
estimated when using AM to define a major merger. 
Since the radius and mass of a subhalo can fluctuate 
due to the temporary incorrect assignment of particles 
to a halo, by using the AM/M > 0.2 criterion more 
than half of the subhalos at z = that were more mas- 
sive than 1O 115 /i _1 M0 at infall, get a merger assigned 
while they were subhalos. However, no secondary pro- 
genitor is found for these subhalos. To overcome this 
problem of incorrect major-merger counting associated 
to the uncertain mass determination of subhalos, we use 
an alternative criterion for subhalos: a major merger is 
counted when M2 mm /M > 0.2, where M^mm is the mass 
of the second most massive progenitor. Note that to 
meet this condition the second most massive progeni- 
tor should have more than 800 particles for halos with 
M > 1O 115 /i _1 M0. Using this criterion, now only 5% 
of all subhalos that survive until z = suffered a major 
merger once they enter in a bigger halo. 

With these definitions we find that halos in high- 
density environments suffered on average more major 
mergers than their counterparts in low-density environ- 
ments. The average number of major mergers, counted 
since the mass of the halo is Mi(a) > 0.05Mi^ max , and 
the standard deviations are 4.8 ±1.4 and 3.9 ±1.3 for 
high- and low- density environments, respectively. The 
difference in the averages is highly significant according 
to Student's t— test. This difference is established mainly 
at z > 3. For all the sample, we measure a mean of 
4.3 ±1.4 major mergers p er halo, similar to the values 
reported in lLi et all l)2005l) . As these authors showed, the 
major merging statistics does not depend significantly on 
halo mass. 

From the pure number of major mergers it is difficult 
to estimate the effect mergers had on a halo, since early 
mergers only contribute a small fraction to the final halo 
mass. Therefore we also measured for every halo the 
fraction of mass accreted in major merger events. For 
the sample including subhalos and in the mass range 
11.5 < logM /h -1 M @ < 12.5, Fig. shows averages of 
the major-merger mass fraction, / mm = M mm /M, vs 
a for all halos (solid line), and for halos in the high- 
(S4 > 5, dashed line) and low- (64 < 0, dotted line) den- 
sity environments. At a given epoch a, M and M mm are 
the current halo mass and the mass assembled in ma- 
jor mergers until a, respectively. The encapsulated panel 
shows the corresponding distribution of / mm at z = Q. 
The denser is the environment, the higher and broader 
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Fig. 7. — Distribution of the ratio of the mass at the last major 
merger Mj mm to the maximum mass M maa; of the halo as a proxy 
for the morphology of the galaxy living in the halo. 
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Fig. 6. — Evolution of the average fraction / mm of mass accreted 
in mergers (AAf/M > 0.2) is shown for all halos in the mass range 
11.5 < logMo//i _1 M < 12.5 as solid line. Halos in high and 
low density environments are shown by dashed and dotted lines, 
respectively. The evolution of / mm for isolated halos is indicated 
by thin lines. The encapsulated panel shows the distributions of 
/mm at z = for all halos. The averages and scatters of /mm at 
z = are given in the main panel. 



distributed is the fraction of mass assembled in major 
mergers. The fraction / mm decreases with a, showing 
that mass accretion becomes more and more dominant 
over major mergers with time. This is particularly evi- 
dent for isolated halos (thin lines in Fig. EJ) . For the high 
density environment / mm remains constant since z ~ 1 
because a significant fraction of halos start to become 
subhalos for which the mass is kept constant. 

Following our results, in dense environments we expect 
more galaxies to assemble early through violent merg- 
ers. However, even in the densest environments, the 
present-day halo mass fraction aggregated through ma- 
jor/medium mergers is small (Fig. EJ). It is commonly 
assumed that the major merger of halos implies the ma- 
jor merger of their corresponding galaxies with a spheroid 
as the outcome. Further mass accretion produces again 
a growing disk. Then we may calculate the ratio of 
the mass when the last halo major merger happened to 
the maximum halo mass, Mi mm /M max , in order to get a 
rough estimate of the spheroid-to-total mass ratio of the 
galaxy. Figure \7\ shows the normalized distribution of 
the Mi, nm /M max ratio for all the halos in the mass range 
11.5 < logMo/Zi^Mg < 12.5 (solid line) at z = 0. Also 
shown are the corresponding normalized distributions for 
the subsamples of halos in the high- and low-density en- 
vironments. For high density environments the distribu- 
tion is shifted to higher Mi mm /M max values. However, 
in all environments almost all halos still accrete at least 
25% of their mass after their last major merger. 

4. DISCUSSION 



We have found that the mass assembly histories of 
galaxy-size CDM halos show marked environmental de- 
pendences, though with large scatters. This kind of 
findings requires numerical simulations due to the com- 
plexity that introduces the local and global environment. 
Popular approaches used for modeling galaxy formation 
and evolution, as the EPS theory and the halo model 
of galaxy clustering, assume that the halo evolution and 
properties depend only on mass but not on environment 
or formation times. The results of the present paper as 
well as of recent wor ks llCao et alJ 120051 IHarker et all 
l2006t IWechslcr et al.ll2005|) should be taken into account 
in these semi-analytical ap proaches (|Sheth fc TormerTI 
l200llAbbas fc ShetHl2T)05f 

In the hierarchical clustering scenario, the evolution 
of galaxies is tightly related to the assembly history of 
their halos. In this section we will discuss and attempt 
to quantify to what extent properties of the observed 
galaxy population are established already at the level of 
CDM halos and how the trends with environment found 
for halos compare to those observed for galaxies. 

Over the last years, a number of groups used 
the technique of grafting semi-analytic models 
of galaxy evolution on to CDM halo merging 
trees constructed from lar ge-volume cosm ol ogical 
N-body simu lations (e.g., iKauffmann et al.l Il999|: 
Diaferio et alJ 120011 iSpringel et alJ 120011: iHellv et alJ 
20031 lHatton et al.l20Q3HDe Lucia. Kauffmann & White! 
2004bt lSpringel et al. ll2005UKang et al.ll2005|) . In these 
works, where the properties, distributions and evolution 
of different galaxy populations were predicted, the 
spatial distribution of galaxies is explicitely taken into 
account. However, in most of these works the halo 
assembly process could only be followed with accuracy 
for the most massive galaxy-sized halos and, in some 
cases, subhalos were not included in the analysis at all. 
On the other hand, it was not quantified explicitely 
to what extent the obtained galaxy properties and 
distributions are the result of the CDM halo assembly 
history. 
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To tackle this question, in the following we will con- 
sider a simple model, in which galaxy SFRs are propor- 
tional to halo aggregation rates and bulge-to-total mass 
ratios (morphologies) are determined by the Mi mm /M max 
ratios of halos. Although we are well aware that this is 
not a realistic model of galaxy formation, the idea is to 
isolate the environmental effects of CDM and to be able 
to compare to observations. 

4.1. The distributions of halo properties 

The recently assembled large galaxy redshift surveys 
such as the Two-degree Field Galaxy Redshift Survey 
(2dFGRS) and the Sloan Digital Sky Survey (SDSS) have 
revealed that s pecific SFRs and colors are not evenly dis- 
tributed (e.g.. Istrateva et al.ll2001t iBlanton et al.ll2003t 
Hogg et al. 1120041 iBrinchmann et al.ll2004t iBalogh et all 
20041 iTanaka et al. I I2004t IKauffmann et al. I 120041: 
WeimnE^iiiet^lTT27)0^. and more references therein). 
Galaxies are clearly separated into two distinct popu- 
lations of red, old and 'passive' early-type galaxies, and 
blue and 'active' star forming late-type galaxies. 

The distribution of the mass aggregation rates and its 
evolution presented in §3.2 seem to lie at the foundation 
of these observed bimodalities of the galaxy population. 
While our simplifying assumption that the specific SFR 
of galaxies is driven mainly by (M / M max )o can in general 
reproduce this bimodality, the aggregation rates at the 
low and high-mass end of the halo population can not be 
reconciled with observations. Here, baryonic processes 
must come into play. The color distribution on the other 
hand is influenced not only by SFRs and ages of the 
stellar population, but also by the metal enrichment, i.e. 
a purely baryonic process. 

One may consider that a fraction of the 'passive' ha- 
los after truncating their MAHs retain for some Gyrs a 
reservoir of hot coronal gas that may still feed the galaxy 
with fresh gas through cooling flows. Therefore, the frac- 
tion of galaxies formed within the 'passive' halos that 
are quiescent should be lo wer than the ~ 40% of 'pas- 
sive' halos found in §3.2. iWeinmann et al.l l)2006j) find 
for the SDSS sample that roughly 31%, 20%, and 48% of 
galaxies belong to their categories of red quiescent, red 
star-forming, and blue actively star-forming galaxies, re- 
spectively. It is interesting to note that these fractions 
could be roughly explained at the level of CDM halo ac- 
tivity if one takes into account that the ~ 20% of red 
star-forming galaxies are formed in some of the 'passive' 
halos and in those with very low (M /M max )o values. 

For the low mass halos which have low aggregation 
rates it is also challenging to expl ain how they can 
host galaxies with high specific SFRs ijBrinchmann et ahl 
2004) . In these halos cooling and star formation must be 
delayed with respect to the dark matter assembly. In 
particular, cold gas may stay in the disk without be- 
ing consumed immediately by star formation. This is 
in general expected, as the SFRs of observ ed galaxies 
are r elated to their gas surface mass density i|Kennicuttl 
1998). For the low mass galaxies, high efficient feed- 
back at low gas disk surface density combined with later 
re-accretion of the gas could result in the high specific 
SFRs at the present epoch. The disk mass surface den- 
sity is indeed p redicted to be lower as the halo mass 
decreases (e.g., iDalcanton. Spergel fc Summers I Il997t 



lAvila-Reese et al. lfl998t IMo. Mao fc Whitdll99"sh . 

We find that the distribution of (M /M max ) is very 
sensitive to environment: 'passive' halos are the major- 
ity in the high density regions, while 'active' halos are 
the majority in the low density regions (see fractions in- 
side Fig. 2}. On the other hand, at least in the mass 
range studied in our simulation, we do not find a strong 
mass dependence in the (M /M max )o distribution (the 
fraction of halos with (M/M max )o < for three mass 
ranges change only by a few percent; see Figs. |3] and 
0J). The main reason of the environmental dependence 
of (M/M max )o is the fact that halos in higher density 
regions become subhalos (which means truncation of the 
mass aggregation process) earlier and much more fre- 
quently than halos in lower density regions. As a result, 
the fractions of 'passive' and 'active' halos in the distri- 
bution of the specific mass aggregation rate changes sig- 
nificantly from low- to high-density environment, while 
the change with mass is only marginal. 

The variation of the (M /M max )o distribution with 
environment indeed resembles the co rresponding varia- 
tions observed for SFR t racers fe.g.. IKauffmann et al. I 
12001 IBalogh et al. II20041 see in particular Croton et al. 
2005, Fig. 2, who use an intermediate-scale environ- 
ment criterion close to the one used here.). However, 
the shapes of the (M /M max )o distributions of halos with 
(M /M max )o > also change with environment. The 
distribution is more peaked in the low-density regions 
than in the high-density ones (Fig. 0J . These differences 
are significant according to a K-S test. The medians of 
M /M max for the halos with positive aggregation rates in 
our two selected environments are 0.02 Gyr -1 and 0.04 
Gyr -1 , respectively. This is at variance with the behav- 
ior that wa s reported for t h e pop ulation of star-forming 
galaxies by IBalogh et al. I {2004) who found no signifi- 
cant differences for these galaxies with environment. 

Observations show a dependence of the bimodality in 
specific SFR or color on luminosity (or stellar mass). It 
is possible that this dependence is due to the correlation 
betw een the galaxy luminosity function and environment 
(e.g.. IHogg et aTll200ll ICroton et al. ll20f)l. The main 
problem for our simplified model arises for the massive 
halos in relatively isolated regions: they aggregate mass 
at a similar or slightly higher rate than the less massive 
isolated halos, in such a way that the sequence of lumi- 
nous red quiescent galaxies in the field is not expected; 
instead an overabundance of luminous blue active galax- 
ies is expected. Observations apparently contradict these 
expectations: there is evidence of a red sequence of lumi- 
nous galaxies in low density environments and luminous 
blue galaxies are rare in any environment ijBalogh et alJ 
2004); it was also found that lower mass field galax- 
ies tend to have higher specific SFRs than higher mass 
field g alaxies, and this tre nd continues to very high red- 
shifts ( Feulncr et alJ2005lL The solution to this apparent 
problem may lie in astrophysical processes such as mass- 
dependent shock heati ng of the halo gas corona coupled 
to AG N feedback fe.g.. lCattaneo et al.l200(j:lBower et alJ 
20061, and references therein). 

4.2. Mass aggregation rate dependence on environment 
and evolution of this dependence 
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Fig. 8. — The figure compares the density dependence of the 
specific mass aggregation rates in the upper two panels at z = 
and 2 = 1 with the dependence of Mi mm /M max ('morphology') on 
density in the lower panels. All halos (isolated + subhalos) in the 
specified mass ranges are included in the plots. 



An interesting question with respect to the environ- 
mental dependences is which halo property shows the 
strongest correlation with environment and if some cor- 
relations are more fundamental than othe rs. Recent in- 
vestigations ijChristlein & Za bludoff 2005) show that the 
SFR-density relation persists, even if stellar mass, mean 
stellar age and morphology are kept fixed. Also color 
and magnitude were shown to have the strongest depen- 
dence on environ ment, compared to su r face-brightness 
and co ncentration (Bla nton et al.ll2f)05|) . iQuintero et al.l 
( 2005) even find that there is no morphology-density re- 
lation at fixed color. To shed some light on the question 
of the relative strength of these effects and a possible 
relation, in Figures 8 and 9 we show the specific mass 
aggregation rates and Mi mm /M max ratios (our proxies 
for the specific SFR and morphology) as a function of 
the density contrast 84 and the cluster-centric radius, 
respectively. 

The overall dependence of the specific mass aggre- 
gation rate, M /M max on the intermediate-scale den- 
sity (64) can be appreciated in the upper panels of Fig. 
03 Solid and dashed lines show the medians of all ha- 
los in the mass ranges 11 < \ogMo/h~ 1 M, 3 < 12 and 
12 < logMo//i _1 M0 < 13, respectively. Dotted lines are 
the 25th and 75th percentiles for the former mass range, 
and the two-dotted-dashed line is the 75th percentile for 
the latter mass range (the 25th percentile falls outside 
the plot). Left panel is at z = and right panel is at 
z = 1. The median of M/M max at z = anti-correlates 
weakly with 64 up to £4 ss 2 — 3; for halos in higher den- 
sity environments, the median of M/M max drops rapidly 
to 0, i.e. the population of 'passive' halos (mostly sub- 
halos) starts to dominate. The differences in these be- 
haviors with mass are small. The more massive halos 



present slightly higher mass aggregation rates. At z = 1, 
there is almost no correlation of the mass aggregation 
rate with environment. Especially halos in high density 
environments are still in their growth phase and no drop 
of accretion rates at these densities is observed. The 
M/M max values at z = 1 are systematically higher by a 
factor of ~ 4 than those found at z = Q. For the most 
massive halos, M/M max even increases slightly with 84. 
The mass aggregation of the latter halos at these epochs 
happens not only by accretion but also, in a significant 
fraction, by major mergers, which are common in a dense 
region that still did not virialize. 

Analysis of the SDSS and 2dFGRS surveys showed that 
SFR tracers, as the Ha line equivalent width, W(Ha), 
correlate with local density of galaxie s or with the 
cluster-centric radius llLewis et al. Il2002t Comez et al. I 
120031 iBaloeh et a!~l 120041: iTanaka et al. 1 12004|) . The 
qualitative features of these correlations are similar to 
those presented here for the mass aggregation rate of 
CDM halos at z = 0. Two populations of galaxies, those 
with significant, on-going SF, and those without SF, are 
revealed. The distributions of I-F(Ha) or the inferred spe- 
cific SFR present an abrupt change at some characteristic 
density p c or cluster-centric radius r c : at densities higher 
than p c (or radii smaller than r c ) there is a near-total 
lack of star-forming galaxies, while at densities smaller 
than p c (radii larger than r c ), the SF activity tends to 
increase. These features are well explained if the current 
SFR of galaxies is associated to the mass infall rate of 
their halos. 

Recent studies show that the galaxy SFR and 
color bimodality is present at redshi fts as high as 



1, though with significant changes llBe l ^t^dJ_|2Qf^ 



Nuiiten et al] 120051 Giallongo et al1l200 
IF 



2005: lCooper et al. 1120 06: Cuc ciati et al. 



.Faber et al. 1 
2006J). Obser- 



vations show that the density of the galaxy red sequence 
roughly duplicates since z ~ 1 to z = with a slight red- 
dening of the peak of the distribution ijBell et al.ll2004t 
iNuiiten et alJl2005l) . From our simulation, we find that 
the fraction of 'passive' halos with M/M max < in- 
creases by factor of ~ 3 since z ~ 1 to z ~ 0. 

To further contrast 'active' and 'passive' halos, in the 
upper panels of Fig. Elwe show the fractions of 'active' 
(M/M max > 0, dotted line) and 'passive' (M/M max < 0, 
solid line) halos in our simulation at z — (left panel) 
and z = 1 (right panel) as a function of cluster-centric 
radius. The use of cluster-centric radius allows to probe 
also the radial dependence inside the group/cluster halo 
where ^4 remains constant. At z = 'passive' halos 
dominate completely inside the virial radius of collapsed 
structures, R V i r . The fraction of 'passive' halos decreases 
beyond R„ r while the fraction of 'active' halos increases. 
In regions as far as ~ 3R V i r , the fraction of 'active' ha- 
los already dominates over that one of 'passive' halos. 
At z = 1 (i) the dependences of the fractions of both 
halo populations on environment become much flatter 
than at z = 0, and (ii) 'active' halos dominate over 
'passive' ones even inside the v irial radii of co l lapsed 
structures. In a very recent work. lCucciati et al. I l)2006fl 
have reported a qualitatively similar evolution for the 
color-density relation by using a sample of thousands of 
galaxies from the VIMOS-VLT Deep Survey. We remark 
that the density criterion used by these authors traces 
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Fig. 9. — The upper two panels show the fractions of 'passive' 
(no accretion) and 'active' (accreting) halos at redshifts and 1 
as a function of cluster— centric radius. The Mi mm /M max ratio is 
symbolized with j3. The lower two panels show fractions of merger 
dominated halos ({3 > 0.5) as solid lines, of accretion dominated 
halos (/3 < 0.3) as dotted lines and of intermediate halos (0.3 < /3 < 
0.5) as dashed lines. Requiring that accretion dominated halos had 
no major merger since z = 1 for halos at z = or in the last 2 Gyrs 
for halos at z = 1 leads to the thin dotted lines. The additional 
condition that the last major merger happened more than 2 Gyrs 
ago for the merger dominated halos leads to the thin solid lines. 



the intermediate-scale environment (~8/i ^-Mpc) rather 
than the local one. 

4.3. The morphology-density relation and its evolution 

The Mi mm /M max ratio distributions for all (isolated + 
sub-) halos in the mass range 11.5 < \ogM /h~ 1 M ( 7 ) < 
12.5, and for the subsamples of halos in the high- 
and low-density environments were presented in Fig. 
[7| The Mi mm /M max ratios are systematically larger 
for higher environmental densities, similar to the ob- 
servational trends for the galaxy bulge-to-total lumi- 
nosity ratio. However, the fraction of halos with large 
Mi mm /M max (thought to be associated to elliptical 
galaxies) is very small, even in the high-density (c l uster) 
environment. This suggests (e.g.. iMaller et alH 12005(1 
that astrophysical processes able to prevent further gas 
accretion onto big spheroids are necessary to explain 
the observed fractions of elliptical and SO galaxies (see 
above) . 

The bottom panels of Fig. |H1 show the median and 
quartiles of the Mi mm /M max distribution as a function of 
64 at z — (left) and z — 1 (right). The line code is as in 
the top panels. At z — 0, there is a weak but continuous 
increase of the M lmm /M max ratio with <5 4 . The correla- 
tion is nearly independent of mass (if any, massive halos 
have slightly larger Mi mm /M max ratios than the less mas- 
sive halos) and it has roughly the same slope for different 
percentiles. Unlike the M/M max distribution, there is 
no sharp change in the values of Mi mm /M max when pass- 



ing from high- to low-density regions. In other words, 
the halo Mi mm /M max -density relation extends smoothly 
from high to low densities. In the case of observed galax- 
ies, the bulge-to-total mass ratio (morphology) indeed 
changes with environment mo re smoothly than the SFR 
l)Christlein fc Z abludofi l2005[) . However, the morphol- 
ogy (mostly of faint galaxies) changes more abruptly at 
a certain density (typic al of clusters outskir ts) than our 
halo Mi mm /M max ratio (Tanaka et al. 112004ft . suggesting 
that baryonic processes should play also a role here. 

At z = 1, the dependence of the current Mi mm /M max - 
ratio on intermediate-scale density is almost absent (see 
lower right panel of Fig. |HJ. This plot can be inter- 
preted as follows: in dense environments (protoclusters), 
the Mi, nm /M max ratio is practically established at these 
early epochs, while in less dense environments, the ha- 
los still grow by mass accretion so that the Mi mm /M max 
ratios decrease with time. The galaxies inside these ha- 
los will grow still blue disks. However, some baryonic 
mechanisms that are effective in the field or poor groups 
can reduce the decrease of the bulge-to-total mass ratio 
(and SFR) a s the density is lowe r. These processes are 
for instance ijTanaka et al. ll2004ft the low velocity inter- 
actions between galaxies, able to trigger SF and consume 
the gas, and strangulation (halo-gas stripping). 

In the bottom panels of Fig. Elthe fractions of merger 
dominated-halos (Mi mm /M max > 0.5), intermediate ha- 
los (0.3 < Mi mm /M max < 0.5), and accretion-dominated 
halos (Mi mm /M max < 0.3) are shown as a function of 
the cluster-centric-radius for redshifts z — (left panel) 
and z = 1 (right panel). The fraction of halos with high 
Mi mm /M max ratios decreases continuously from ~ 0.4 in 
the center of virialized host halos to ~ 0.2 at 4 — 5R V i r . 
The fraction of halos with low Mi mm /M max ratios is small 
in the centers of virialized host halos (~ 0.25), but it con- 
tinuously increases from the periphery to regions of low 
density where it attains a value of ~ 0.5. The fractions 
of halos with low and high Mi mm /M max ratios become 
comparable at ~ lR V i r . Intermediate halos are more 
abundant inside virialized host halos than in low-denstiy 
regions. 

The observed morpholog y-density (e.g.. iDressler 
1980f iPostman fc Geller I 119841: IDressler et aT 
l!H IQoto et a,l. I I20Q1 



centric radius 



iDommguez et al. 112001 



e.g. 



and morphologry- cl uster- 



Whi tmore fc GallagheH Il99ll 
G^to^t*ldT]2003(l relations for 



galaxies are in qualitative agreement with the trends ob- 
tained here for halos. The fraction of intermediate-type 
galaxies (mostly SO's) increases from the outskirts of 
clusters to intermediate cluster radii, while the fraction 
of late-type galaxies decreases. In the densest regions 
(cluster centers), the intermediate- and late-type 
galaxy fractions decrease, while the early-type fraction 
increases. In the sparse regions (far from clusters), the 
morphology-density relation flattens, i.e. the fractions of 
different galaxy types only slightly change with density 
or cluster-centric radius, and the fraction of late-type 
galaxies dominates. Probably, the main difference with 
halos is that the fractions of merger- and accretion- 
dominated halos continue depending on cluster-centric 
radius (or on density) in these low-density environments 
(see Fig. EJ. However, there are several baryonic 
processes that could work in the correct direction to 
"flatten" the morphology-density relation (see above). 



DEPENDENCES OF CDM HALO ASSEMBLY ON ENVIRONMENT 



11 



At z = 1 the fraction of halos with high Mi mm /M max 
values in the inner parts of clusters is higher (~ 0.6) 
than at z — (~ 0.4). The existing groups/clusters at 
z = 1 just formed and contain still a large fraction of 
halos that suffered a major merger very recently (large 
M] mm /M max ratios). At z = a large part of those halos 
has merged with the parent halo and is replaced with 
new subhalos that had more time to aggregate mass by 
smooth accretion before they became subhalos, resulting 
in a higher fraction of halos with smaller Mi, nm /M max 
ratios than at higher redshift. In regions outside virial- 
ized groups/clusters, most of the halos continue accreting 
mass until the present epoch, so that their Mi mm /M max 
ratio decreases more and more, increasing therefore dras- 
tically the fraction of halos with small Mi mm /M max ratios 
(compare also Fig. 6). 

When attempting to connect our results with obser- 
vations of galaxy populations, we should have in mind 
that the halo last major merger may have happened 
very recently. In this case a high Mi mm /M max value 
does not imply a residing early-type galaxy, but rather 
a star-bursting one. If we impose the extra condition 
that the last major merger happened more than 2 Gyrs 
ago, then the fraction of merger-dominated halos de- 
creases, mainly at z = 1 (thin solid lines in lower panels 
of Fig. EJ. Nevertheless, this fraction is still larger at 
z = 1 than at z = 0, contr ary to the obsery ational infer- 
ences l)Smith et al. 1120051 IPostman et al. Il2005|) . Since 
the "morphological" fractions reported here for halos at 
z = are comparable to those observed for galaxies, the 
main disagreement actually occurs at z = 1. 

On one hand, at high redshifts and in dense environ- 
ments, most of the halos are still actively growing, with a 
significant mass fraction aggregated in mergers. There- 
fore, a part of the high-density halo population with 
Mi mm /M max > 0.5 at z — 1 is expected to harbour actu- 
ally blue, star-forming galaxies. On the other hand, the 
strong increase of the fraction of accretion-dominated 
halos since z — 1 to z = 0, especially at large cluster- 
centric radii, contrasts with the correspon ding decrease 
in th e observed late-type g alaxy fraction l)Smith et al. I 
120051 IPostman et al. 1 120051. As discussed above, bary- 
onic processes in the field and poor groups may work to 
decrease the present-day fraction of late-type galaxies 
in the low-density accretion-dominated halos. However, 
it is not easy to explain the observed high fraction of 
late-type galaxies in low-density environments at high 
redshifts as compared to the low fraction of accretion- 
dominated halos in these environments. Selection ef- 
fects could be affecting the observational inferences of the 
morphological mix. Detailed galaxy modeling inside the 
evolving ACDM halos is necessary to understand better 
the evolution of the morphology -cluster-centric radius 
and -density relations. 

Finally, by comparing the strenght of the environ- 
mental dependence of the aggregation rate and the 
Mi mm /M max ratio presented in Figures 8 and 9, we con- 
clude that Mi mm /M max shows a much weaker dependence 
on environment than the aggregation rate, M / M max . 
It could be, however, that the difference between the 
strength of the observed morphology-density relation 
and our Mi mm /M max -density relation is caused by pro- 
cesses influencing both, SFRs and morphology, at the 



same time. 



5. CONCLUSIONS 



In previous works it was shown that the 
mass function jLgmggi^^Kauffmann 



nal properties (Avila-Recsc et al. 2005) of ACDM halos 



resent-day 
and inter- 



change systematically with environment. In the present 
paper we investigated environmental dependences of the 
mass assembly process of ACDM halos. Our results 
suggest that a significant part of the relations observed 
between galaxy properties and environment can be ex- 
plained by the environmental dependences of the halo 
assembling process. Our main findings are as follows: 

• The MAHs of the ACDM halos change with 
intermediate-scale environment: the lower the density 
contrast ^4, the later on average a given fraction of the 
maximum halo mass is assembled. The averages of the 
specific present-day mass aggregation rate, (M /M max ) , 
are 4-5 times higher for halos in low-density environ- 
ments (84 < 0) than for halos in high-density environ- 
ments (84 > 5). These differences are smaller at higher 
redshifts, at z ~ 1 disappear, and for z > 1 the trends 
are reversed. The average MAHs and mass aggregation 
rates also depend on halo mass, but to a much lesser 
degree than on environment. 

• The distribution of (M /M max )o has two parts: 
w 40% of all the halos in the mass range 11.5 < 
log M o /ft. _1 M < 12.5 do not aggregate mass, i.e 
(M/M max ) < 0, while the rest has a distribution of 
(M/M max ) peaked at - 0.03 - 0.04 Gyr" 1 . The dis- 
tribution changes little with mass but dramatically with 
environment: the fraction of halos with (M/M max )o < 
increases to w 81% and reduces to « 22% for the high- 
and low-density environments, respectively. The distri- 
bution of the halo formation time, also changes with en- 
vironment: for instance, halos of ~ 10 12 h~ 1 MQ assemble 
half of their maximum mass approximately 1.5 Gyrs ear- 
lier in the high-density regions than in the low-density 
ones. 

• Using only isolated halos, the dependences of halo 
MAHs on environment are weak, showing that the 
intermediate-scale density around halos affects only par- 
tially their mass assembling process. The strongest ef- 
fects of environment appear when subhalos are included. 
Although the effect due to subhalos is local, the fraction 
of halos with masses above 10 11 /i _1 M s that become sub- 
halos increases as the environment gets denser. The two 
main effects that a halo suffers when it becomes subhalo 
of a larger host virialized structure are: (i) its MAH is 
suddenly truncated (the mass growth is even reversed to 
mass loss due to tidal stripping), and (ii) the major merg- 
ing probability drops drastically due to the high velocity 
dispersions inside the virialized host halo. 

• Present-day halos in high density environments suf- 
fered more major mergers on average and assembled a 
larger fraction of their mass in major mergers than halos 
in low-density regions. Halos in dense regions become 
subhalos at z < 1 much more frequently than those in 
low-density regions. For these halos the (high) major 
merging mass fractions and (large) last major merger 
mass-to-current mass ratios, Mi mm /M max , become de- 
fined at the epoch of their infall. On the other hand, ha- 
los in low-density environments continue growing mainly 
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by mass accretion, thus, their major merging mass frac- 
tion and Mi mm /M m ax ratio continuously decrease. 

We have shown that the distribution of (M/M max )o 
and the strong dependence this distribution presents 
with environment could explain partly the observed bi- 
modal distributions of specific SFR and color for galax- 
ies and their systematical dependences on environment. 
The main difficulty for this direct halo-galaxy connection 
would be to explain the existence of the population of red 
luminous galaxies in low-density environments, and the 
high fraction of blue star forming galaxies among low 
mass galaxies. Astrophysical processes related to bary- 
onic matter should be invoked here. 

The morphology-density relation can be partly ex- 
plained if the halo Mi mm /M max ratio is used as measure 
of the bulge-to-total mass ratio (morphology). However, 
even close to the cluster center, the fraction of halos with 
Mimm/Mmax as high as to harbour elliptical galaxies is 
too small at z = 0. This suggests the need of astro- 
physical processes able to avoid gas accretion onto the 
early formed spheroids. The main trends of the observed 
morphology- density and cluster-centric radius relations 



at z = agree with the trends of the halo Mi mm /M max 
ratio with environment reported here. 

Future comparisons of observations with models that 
include the astrophysical processes of galaxy formation 
in highly resolved evolving ACDM halos will tell us 
whether or not the ACDM paradigm agrees with obser- 
vations, and what are the key ingredients at the basis 
of the morphology-, color-, and star formation-density 
correlations. Our results suggest that the halo assem- 
bly process and its dependence on environment should 
be one main ingredient. 
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